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In this letter the paper [R. Aquilano, M. Castagnino, Mod. Phys. Lett A, 11, 755 (1996)] is 
improved by considering that the main source of entropy production are the stars photospheres. 



I. INTRODUCTION. 

In papers jjj Q and Q two of us reported a rough coincidence between the time where the minimum of the entropy 
gap A5* = Sact — Smax |§ , takes place and the time where all the stars will exhaust their fuel. The rough nature 
of this result was explained saying that the calculation was done using a model extremely naive and simplified: a 
homogeneous universe In fact in the real universe nuclear reactions (that were considered as the main source 
of entropy) take place within the stars, that can only be properly considered in an inhomogcneous geometry. So 
we expend some time trying to generalize our model to these type of geometries and computing the entropy of the 
galaxies and stars in the condensation period, after decoupling time, in order to make the coincidence precise ^. It 
was soon clear to us that the latter entropy can be neglected with respect to the amount of entropy produced within 
the stars and therefore the mechanism that produced the inhomogeneity is not so important for our problem. Stars 
are formed when a sufficiently large mass of interstellar gas is compressed into a small enough volume, so its force of 
self-gravitation becomes sufficiently great to cause gravitational collapse. The instability makes the cloud to reduce 
its size quickly, the temperature rises and pressure forces begin to restrict the collapse until hydrostatic equilibrium is 
obtained. In this process, gravitational energy is converted in kinetic energy and radiation M]. In the classical model 
of star formation the time of contraction is the Kelvin-Helmholtz time || : 

GM 2 

where M is the mass, R is the radius and L is the luminosity of the star. For a star like the sun this time is 3 10 7 
years || . In more recent models of star formation, the relevant time scale is the much shorter free-fall time || : 

t ff = (pG)-^ = 7 -lO 5 ^) years (2) 
where p and n are the cloud mass density and the cloud number density respectively. 

In any case, it takes less than 1% of the lifetime (= Gyr) to the star to form. In the contraction to make the 
star, the energy radiated is a part of the gravitational energy and the nuclear reactions are not important. The 
binding gravitational energy is negligible compared with the energy liberated through nuclear reactions when the star 
is already formed Q . As a consequence we can neglect the entropy production in the phase of contraction that makes 
the star. 

So, in order to improve our rough result, we were forced to change the approach. 



1 It is a good occasion to correct some errata in the abstract of paper |Tj . 
-In line 2: it reads "quantitative", it must read "qualitative". 
-In line 3: it reads "entropy", it must read "entropy gap". 
-In line 4: it reads "qualitative", it must read "quantitative". 
Besides that the higher order terms in eq. (3) of ref. Jl) were neglected; perhaps they may be important for finite times. 
3 In fact, inhomogeneity must be considered, but in a different way, as we will see below. 
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Let us review the main formulae. The entropy gap was the conditional entropy of p(t), the state of the universe at 
time t with respect to the equilibrium state at that time p*(t) : 



AS = -tr[plog(p- 1 p)] (3) 

where 

p(t) = p„(t) + Pl exp(-Tt) + Ol^t)- 1 } (4) 

and 

p,M = ZT- 3 -J—- (5) 

exp( T - 1) 

pi was a phenomenological coefficient constant in time, the time variation of the main irreversible process was 
exp(— jt) so 7 _1 was the characteristic time of this process, T was the temperature of the universe, T ~ a -1 where a 
is the radius or scale of the universe, and Z is a normalization constant. 

Then we could approximate AS as 

AS = -CT 3 cxp(- 7 i)exp(^) (6) 

where uj\ was the characteristic energy where pi is peaked. 
The time where the minimum of AS* is located was: 

We selected numerical values for the parameters: loi = Tj^ Rl the temperature of the nuclear reactions within the 
stars (that was used in paper JIJ as the main source of entropy), t^R = 7 _1 the characteristic time of these nuclear 
reactions, to the age of the universe, and To, the cosmic micro- wave background temperature. 

Tvr and t were chosen between the following values [|| : 

T NR = 10 6 to 10 8 K (8) 
t NR = 10 6 to 10 9 years 

while for to and To we can take: 

t Q = 1.5 x 10 10 years (9) 
To = 3K 

In order to obtain a reasonable result we chose (with no explanation) the lower bounds for Tnr and t^R and for t cr 
we obtained : 

t cr d 10% - I0 u years (10) 

concluding that the order of magnitude of t cr was a realistic one. In fact, 10 4 io w 1-5 x 10 14 years after the big-bang 
all the stars will exhaust their fuel |Qq so it is reasonable that this time would be of the same order than the one 
where the entropy gap stops its decreasing and begins to grow Q]. But the choice of the lower bound in eqs. (^) was 
not explained and only the inhomogeneity was argued as above. 



'After the publication of paper in 1996, new data about this time appear in 
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II. THE PHOTOSPHERE AS THE MAIN UNSTABLE SYSTEM. 



Now we have reconsidered the problem and we conclude that, even if nuclear reactions within the stars are the main 
source of entropy, the parameters Tnr and tNR are not the good ones to define the behavior of the term exp(— ft)pi 
of equation (||)((3) of paper since they do not correspond to the main unstable system that we must consider [| 
In fact the main production of entropy in a star is not located in its core, where the temperature is almost constant 
(and equal to Tjvij), but in the photosphere where the star radiates. The energy radiated from the surface of the star 
is produced in the interior by fusion of light nuclei into heavier nuclei. Most stellar structures are essentially static, so 
the power radiated is supplied at the same rate by these exothermic nuclear reactions that take place near the center 
of the star Once the star is formed, it settles into a termally stable state where all the nuclear energy is radiated 
at the surface and the rate of internal entropy change is extremely low ||. 

We can decompose the whole star in two branch systems ||, as explained in paper Q| (or in section VII of paper Q), 
where a chain of branch systems was introduced. We have two branch systems to study: the core and the photosphere. 
The core gives energy to the photosphere and in turn the photosphere diffuses this energy to the surroundings of the 
star, namely in the bath of microwave radiation at temperature To .In this way, we have two sources of entropy 
production: the radiation of energy at the surface of the star and the change of composition inside the star (as 
time passes we have more helium and less hydrogen). Since the core of a star is near thermodynamic equilibrium, 
we neglect the second and we concentrate on the first: the radiation from the surface of the star (related with the 
difference between the star and the background temperatures) . So the temperature of the photosphere and not the 
one of the core must be introduced in our formula. This is also the case for the lifetime. We must take the lifetime 
of the photosphere not the one of the nuclear reactions. Thus it is better to consider the photosphere as the unstable 
system that defines the term exp(— jt)pi of equation (^). So we must change Tnr and t^R by Tp, the temperature 
of the photosphere and ts the characteristic lifetime of the star respectively. Then we must change eq. (0) to: 

/2T P t s \ 3 , . 

Considering that the mean mass of stars is 0.64M© with surface temperature 4.6 • 10 3 K and lifetime 3.8 • 10 10 years 
(see appendix) we obtain: 

t cr = I0 10 t = 10 20 years (12) 

but now the computation was not done using an arbitrary choice of the lower bound in some data, but using the 
first meaningful figure in all the data []. 



III. THE STELLIFEROUS AND THE DEGENERATED ERAS. 

Let us now compare this result with those of paper |HJ|, where the future history of the universe is analyzed. 
t cr = 10 10 io = 10 20 years is placed after the end of the " stelliferous era" (10 6 < t < I0 14 years), where most of the 
energy generated in the universe arises from nuclear processes in the conventional star evolution, and at the beginning 
of the "degenerated era" (10 15 < t < I0 37 years), where most of the (baryonic) mass of the universe is locked up in 
degenerated stellar objects: brown dwarfs, white dwarfs, and neutron stars. In this era energy is generated through 
proton decay and particle annihilation. 

So t cr = 10 20 years is not a bad place for the minimum of the entropy gap, since it is bigger than the end of 
conventional star formation (10 1 years), it is also bigger than the typical time of star formation via brown dwarfs 
collision (10 16 years), and it is of the order of the time corresponding to stellar evaporation from galaxies (I0 19 years). 
Namely a time where we can consider that the main mechanisms of formation of stars are ended while the beginning 



5 In fact, to take the parameters Tnr and ijvjj corresponds to take an homogenous gas model that fills the universe where, 
nuclear reactions take place in this homogenous gas, while in real universe these reactions take place within the star in a quite 
inhomogeneous scenario. 

6 In paper JhJ, using the photosphere data, but just orders of magnitude, we have again obtained (p"ci|). 

7 We have obtained a larger result than (^) , but remember that the latter was obtained just by taking the lower bound in eqs. 
(^). If we would choose the mean values in these equations we would obtain a larger result than the one of eq. dl^), without 
any significance. 
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of the evaporation process has just started. This is the best place for t cr , that must be the frontier between the 
growing order period (formation of structures with decaying entropy |J) and the diminishing order period (decaying 
of the structures with growing entropy). Tacking into account the great uncertainty of all cosmological data we can 
say that t cr is located in the edge between the stelliferous normal era, where we are living and which is dominated 
by the formation of structures, and the future degenerated era, full of, by now strange, objects, where the growing of 
disorder will begin. So even if eq. (jll]) is the result of a very simple model it gives a very reasonable value. 



IV. CONCLUSION. 



By choosing a more realistic model as the main source of entropy production, the photosphere of the stars, we have 
obtained a reasonable value for the time where the minimum of the entropy gap is reached. This is the frontier of 
the formation and destruction of structure periods and therefore one of the most important moments of the future 
history of the Universe. 
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VI. APPENDIX 



The stars have masses in the interval [[12| , |T^] : 

O.K m = — - < 100 

M " 

where M» is the mass of the star and M© is the mass of the Sun. 

The mean mass of stars can be calculated using the initial mass function (IMF) £(m) defined by jll 

dN = £(m) • dm 



where dN is the number of stars with masses between m and m + dm 
We will use as IMF [0: 



dN 
dQxim] 



= ip(hi m) — exp 



A - 



2<a>< 



■In" — 



(13) 



(14) 



(15) 



where < a >= 1.57 , m c = 0.1 and A is a constant that sets the overall normalization of the distribution (which 
is not important because it cancels in the calculation of the mean mass of stars). This form of the IMF is consistent 
with observations (see pi). As 



a s m \ d(lnm) 1 1 

f(m) = ip[nim) ■ — : = — ip(mm) = — exp 

dm m m 



A — 



1 



2 < (7 > 2 



In 2 



(16) 



then the mean mass of stars is: 



< 771 > = 



100 

J TTi • • dm 

OA 

100 

J £(m) • dm 
o.i 



= 0.64 



(17) 



Let us calculate the effective surface temperature of a star with mass m — 0.64. The Luminosity of a star is given 
by the Stefan- Boltzmann law: 



L = A-kR oT2 



(18) 
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where R is the radius of the star, a is the Stefan- Boltzmann constant and T p is the effective temperature of the 
surface of the star. Then 



L ( R \ 2 / T„ x ; 



Lq \ Rq J \T ( 

where Lq , Rq ,Tq are, respectively, the luminosity, the radius and the surface temperature of the Sun. 
The radius and luminosity for stars in the main sequence, as a function of mass are (see flq]): 



R _ , / M * x ' 



m? (20) 



(19) 



(21) 



l_ s (M m y 

Lq ~ [Mq) 

with f3 = 0.6 for stars of low mass and r\ = 3.2 .So we have 

Tp-T ^) = Tq • m ' 5 (22) 

Using that T© = 5780/^ and m = 0.64 we obtain 

T P = 4.6 • lO 3 ^ (23) 
The lifetime of a star can be calculated by the equation(see 0] ) : 

( M \ ~ a 

t s = 10 10 ( j^- \ years = 10 10 m- a years (24) 

where a = 3 — 4 for stars of low mass. Taking a = 3 we have: 

t s 3.8 • 10 10 years (25) 
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